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Abstract 

We consider interacting one-dimensional, spinless Fermi gases, whose low- 
energy properties are described by Luttinger liquid theory. We perform a 
systematic, in-depth analysis of the relation between the macroscopic, phe¬ 
nomenological parameters of Luttinger liquid effective field theory, and the 
microscopic interactions of the Fermi gas. I11 particular, we begin by ex¬ 
plaining how to model effective interactions in one dimension, which we then 
apply to the main forward scattering channel - the interbranch collisions - 
common to these systems. We renormalise the corresponding interbranch 
phenomenological constants in favour of scattering phase shifts. Interest¬ 
ingly, our renormalisation procedure shows (i) how Luttinger’s model arises 
in a completely natural way - and not as a convenient approximation - from 
Tomonaga’s model, and (ii) the reasons behind the interbranch coupling con¬ 
stant remaining unrenormalised in Luttinger’s model. We then consider the 
so-called intrabranch processes, whose phenomenological coupling constant 
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is known to be fixed by charge conservation, but whose microscopic origin is 
not well understood. We show that, contrary to general belief and common 
sense, the intrabranch interactions appearing in Luttinger liquid theory do 
not correspond to an intrabranch scattering channel, nor an energy shift due 
to intrabranch interactions, in the microscopic theory. Instead, they are due 
to interbranch processes. We finally apply our results to a particular example 
of an exactly solvable model, namely the fermionic dual to the Lieb-Liniger 
model in the Tonks-Girardeau and super- 1 Tonks-Girardeau regimes. 

Keywords: Luttinger liquids, Renormalisation, Effective field theory, 
One-dimensional systems 


1. Introduction 


Interacting many-body systems in one spatial dimension, once considered 
theoretical playgrounds for the more involved “physical” three-dimensional 
systems, hold currently, and for quite some time, the status of physically rel¬ 
evant theories. Effective reduced-dimensional systems are indeed routinely 
prepared and studied experimentally. Prominent examples include ultracold 
atomic gases 0 , 0, 0 B, 0,0,0, organic conductors [3, 0, 0 and nan¬ 
otubes llj. From the theoretical point of view, one-dimensional systems 
are very appealing, for a number of reasons. Firstly, there are a number of 
exactly solvable models in one spatial dimension. Examples include those 
models solved via the Bethe ansatz technique 0], such as the spinless Bose 
gas with Dirac delta interactions (Lieb-Liniger model) [13j and its spin-1/2 
fermionic counterpart 0. and models with unbroken SuSy, which admit 
Jastrow-product ground state wave functions, such as Sutherland’s model 
with inverse-square interactions and variations thereof [0. More recently, 
strongly-coupled mesoscopic multicomponent systems in arbitrary external 
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potentials have been the subject of intense investigation |16|, Ira 
Moreover, some numerical methods work especially well for one-dimensional 
systems. For instance, the density matrix renormalization group (DMRG) 
|20| is extremely powerful and flexible to use in most one-dimensional many- 
body problems on a lattice. 

The low-energy properties of many of the above systems can be described 


by the universality class of Luttinger liquids |23j, l24j . In brief, the proper¬ 


ties of these systems, such as the collective excitation spectrum, correlation 
functions, or density of states, at energies close to the Fermi energy, exhibit 
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universal behaviour that is described by the exact solution of Tomonaga 


Luttinger’s model 2l|, 22J via bosonization |25|, 126), 127|, l28|, [29|]. For instance, 


the density of states behaves as a power law, where the exponent and the pro¬ 
portionality constant are the only microscopic, non-universal details, around 


the Fermi points 301, l3l|, |32 


Historically, there are two fundamentally different, yet complementary 
approaches for the introduction of Luttinger Liquid phenomenological mod¬ 
els for one-dimensional quantum gases. The first one, called constructive 
bosonisation [33ij], naturally arises for fermionic systems upon linearisation of 
the kinetic energy dispersion, as introduced by Tomonaga [2l[. The second, 
more phenomenological, and also more general approach, is field-theoretical 


bosonisation, introduced by Haldane [34 


In this article, we consider one-dimensional spinless Fermi gases in the 
fermionic Luttinger liquid regime, that is, when their low-energy behaviour 
can be modelled by Luttinger’s model. We will follow the constructive, 
rather than the field-theoretical route, which is most convenient for weakly- 
interacting fermions. In particular, we study the microscopic origin of the 
interaction coupling constants in Luttinger’s model. We first address the 
main forward-scattering channel in the microscopic theory, corresponding to 
the collision of two fermions with incident momenta close to the two op¬ 
posite Fermi points. In the very weak-coupling regime, these interbranch 
processes are well understood, and we show how to renormalise Luttinger’s 
interbranch interaction in favour of scattering data (i.e. the two-body phase 
shift). From our analysis of the interbranch collision channel, we show how 
Luttinger’s model naturally arises from Tomonaga’s model upon renormali- 
sation of the Fermi points, which are to be considered as bare parameters of 
the theory. In this way, we also explain the reason behind the interbranch 
coupling constant remaining ”unrenormalised” in Luttinger’s model, a well 
known fact, as opposed to Tomonaga’s model, where the theory requires 
non-trivial renormalisation. We then consider intrabranch interactions, i.e. 
interactions between particles moving in the same direction. In this case, 
symmetry dictates that the intrabranch coupling constant must be identical 
to the interbranch one [35]. However, the current understanding of the mi¬ 
croscopic origin of the intrabranch interaction in Luttinger’s model is quite 
unsatisfactory. Usually [35;, 36], its origin is said to be linked to the forward 
scattering of fermions near the same Fermi point, i.e. in the same branch. 
Even though this picture sounds reasonable, and easy to grasp physically, it 
has some major issues. Here, we carefully illustrate these issues, and find 
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that the effective mtrabranch interaction in Luttinger’s model corresponds, 
in the microscopic fermionic model, to an energy shift due to mierbranch 
interactions. We finally apply our results in order to obtain the Luttinger’s 


parameter and the speed of excitations in the Tonks-Girardeau regime 37| of 


the Licb-Liniger gas, corresponding to the weak-coupling limit of its fermionic 
dual 
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2. Microscopic system 

We consider a system of N non-relativistic spinless fermions of mass m 
interacting via pairwise potentials W in one spatial dimension (ID). The 
Hamiltonian in the first quantisation has the form 

fr2 N X2 N 

h = ~ 2 ;s£a^+« 

i=l 1 i<j =1 

where X{ is the position of the i-th particle. In the momentum representation, 
the two-body interaction is given by 

/ OO 

dxe iqx W(x), (2) 

-OO 

where we have defined q = k — k'. With this notation, we can write down 
the second-quantised Hamiltonian in the momentum representation, 

H = e ( fc )4 Cfc + ^ V (®) c l+ q c l- q c k'C k - (3) 

k kk'q 

Above, e(k) = h 2 k 2 /2m is the single-particle kinetic energy dispersion, L is 
the system’s size, C& (c£) annihilates (creates) a fermion with momentum k, 
and therefore satisfies canonical anticommutation relations 

{ ( 'k ■ ^k.k/ ■ ( 4 ) 

In Hamiltonian (J3J) we have assumed periodic boimdary conditions, that is, 
k = 2‘Kn/L with n being integer numbers. 

Since we are working with fermions, it is convenient to introduce even 
and odd wave components in the interaction as h = V s + V p . The even 
(odd) wave component, acting only on bosonic (fermionic) wave functions, 
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sometimes called ”s-wave” (”p-wave”), is given by the following symmetric 
(antisymmetric) combination 


V s (k,k') = l[V(k,k') + V(k,-k')], 

(5) 

V p (k,k , ) = ^[V(k,k')-V(k,-k')}. 

(6) 


In the following, we will use V or V p at our discretion, depending on whichever 
is more convenient in each situation. 

3. Model interactions 

In many applications, it is advantageous to construct effective interactions 
that mimic the exact scattering properties of a system at the relevant energy 
scales. This can greatly simplify the calculations while keeping a correct de¬ 
scription of collisional processes, and allow for universal (i.e. independent of 
microscopic details) predictions of many quantities of interest. For instance, 
low-energy scattering is adequately described by means of ’’pionless” effective 
field theory (EFT) [39[], in which the interaction is Taylor-expanded around 
q = k — k' = 0 as 


V(k,k') = g a + g b (k-k') 2 + .... (7) 


The lowest-order contribution to the p-wave potential above is given by 

V p {k, k') « -2 g b kk' = gkk' = V p {4 \k, k'), (8) 


where g is the bare coupling constant that is to be renormalised in favour of 
low-energy scattering data. 

In this article, we shall be interested in collisions of particles with relative 
momenta of 0(fcp), where k F is the Fermi momentum. The EFT interaction 
(HD can be used if the particles have small relative momenta kk'/k F <C 1. If, 
instead, kk' = 0(k F ), the lowest order interaction is given by 


V p (±k F ,±k F )n -[V(0)-V(2k F )\ 


(9) 


The resulting low-energy interaction is sometimes referred to as ’’constant” 
23]. This is a misconception that can lead to contradictions [35]. If we wish 
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Figure 1: Left: depiction of interbranch scattering processes around the Fermi points. The 
red solid line corresponds to the Galilean (parabolic) dispersion, and the blue dashed line 
corresponds to the linearised approximation around the Fermi points. Right: interbranch 
interaction diagram for spinless fermions. 

to use Eq. ([9]) for all momenta, this must satisfy V p (k, —k') = —V p (k, k') (see 
Eq. ()6j)). This gives a model interaction of the form 

V p (k, k') « ^ [E(0) -1/(2 k F )\ sgn(A:)sgn(fc / ) 

= f^sgn^sgn^') = 1/ (2) (/c, k'). (10) 

The collisions due to the above interactions, and its corresponding diagram, 
are depicted in Fig. |T| For the sake of completeness, we write down the 
interaction V p ' in the second-quantised form 

V i 2) = H E S ^~ k - + ^+k +i C Q-k-q C Q-^Q+^ ( n ) 

Q,k,q 

where we have relabelled the momenta in the sum as Q = K/2 = {k + k')/ 2, 
k = (k — k')/ 2 and q = q/2. 

4. Collision-theoretic treatment 

All relevant scattering properties of the system can be extracted from the 
transition matrix (T-matrix) T[z). After separation of centre of mass and 
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relative coordinates, the T-matrix for an odd-wave interaction satisfies the 
following Lippmann-Schwinger equation 


T(z) = V„ + V p G„(z)T(z), 


( 12 ) 


where z is a complex (relative) energy, and Gq(z) = [E — H 0 ]~ l is the non¬ 
interacting Green’s function, with H 0 the non-interacting relative Hamilto¬ 
nian. The scattering states of a Hamiltonian of the form H = H 0 + V 
correspond to an incident wave | k) with non-interacting energy e(k) and a 
scattered wave \^f k )- This is given by 

\^f k ) = lim G 0 (e(fc)+i//)T(e(fc)+ir/)|fc) = G 0 (e{k)+M ± )T{e{k)+ i^)\k). 

r/->0± 

(13) 

Above, the ± signs correspond to waves propagating forward (+) and back¬ 
wards (—) in time, as dictated by time-dependent scattering theory [i(|. The 
non-interacting Green’s function on the energy shell, Go(e(k) + ir /), can be 
written as 


(q\G 0 (e(k) +iO ± )|g) = 2tu 5(q 


V 


e(k) - e{q) 


=F 27rh(e(/c) 


e(?)) 


(14) 

where V denotes Cauchy’s principal value. We now define the principal value 
T-matrix T and the principal value non-interacting Green’s function Qq as 


Qo(e(k)) — — [Go(e(/c) + i0 + ) + Go(e(k) — iO )] , (15) 

T(e(k)) = i [T(e(k) + i0+) + T(e(k) - i0“)] . (16) 

It is easy to see that the principal value T-matrix 0 for real positive energies 
E satisfies the Lippmann-Schwinger equation 


T(E) = V p + V p g 0 (E)T(E). (17) 

The advantage of using the principal value T-matrix is that it avoids the 
use of imaginary - and irrelevant except for bound states - parts, and it 
completely determines the scattering properties of a system of two identical 
particles. 


lr This is also called reaction matrix (l?-matrix) in the literature. 
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4-1. Interbranch scattering: quadratic dispersion 

We now solve the Lippmann-Schwinger equation, Eq. urn for the inter- 

/ o\ 

action Vp J in Eq. (TTU]) . The T-matrix is readily calculated and reads 

T(E; k, k') = g 2 sgn(k)sgn(k'). (18) 


After straightforward but tedious integrations, we find the position-represented 
fermionic wave functions in the relative coordinate x for relative momentum 
k and energy h 2 k 2 /2p (p = m /2 is the reduced mass), 


^k{x) 


1 


mg 2 

nh 2 1 k | 


Ci(|/cx|) 


sin {kx) 


-^jUsi(fcx) cos (kx). 


(19) 


Above, Ci (Si) is the cosine (sine) integral. The only relevant information in 
the above wave function is the asymptotic (x —> ±oo) behaviour. At long 
distances, the wave function in Eq. (fl9|) behaves as 


ipk( x ) oc sin (kx) 


mg 2 
2 h 2 k 


sgn(x) cos (kx). 


( 20 ) 


We compare the above asymptotic wave function with the general asymptotic 
form of a fermionic scattering state in ID, ^A(x) oc sgn(x) sin(/c|x| +9k), and 
we immediately identify the relation between the coupling constant g 2 and 
the scattering phase shift 9 k , as 


tan 9 k 


mg 2 
2 h 2 k' 


( 21 ) 


We still need to decide upon a renormalisation condition, that is, the colli- 
sional energy at which we wish to renormalise the interaction, in which case 
the phase shift at that specific energy is exact. Naturally, we should renor¬ 
malise the interbranch scattering phase shift for particles near opposite Fermi 
points ±/c_f, that is, with relative momentum kp- This yields the following 
value for the bare coupling constant 


92 = 


2 h 2 kp 
m 


tan 9 kp 


—2 hvp tan 9 kF i 


( 22 ) 


where we have defined the Fermi velocity vp = hkp/m. 









4-2. Interbranch scattering: linearised dispersion 

We now consider the approximate linearised spectrum around the Fermi 
momenta ±kp\ which is the starting point for building Luttinger liquid 
theories for fermions. We have included the ® superscript since we would 
like to treat, in general, the Fermi points as bare quantities of the theory. 
This might sound unusual, but notice that (i) upon any approximation all 
coupling constants of the theory may need to be renormalised and in this 
case it is the dispersion itself that has been expanded (and truncated) as a 
power series; and (ii) renormalisation of the Fermi points is always implicitly 
assumed in Luttinger liquid phenomenology j23|, as we shall see below. 

The single-particle energy dispersion is well approximated by 

e{k) fa + hv F (\k\ - kP), (23) 

provided that | k — k^\/k^ <C 1. In Eq. (j25j) we have defined the (bare) 
Fermi energy E^ = h 2 [kp^] 2 /2m. The many-fermion system with the lin¬ 
earised dispersion (1231) corresponds to Tomonaga’s model j4lj, where k e 
(—oo,oo). In order to connect with the phenomenological model later, we 
will call right (left) movers particles with k > 0 (k < 0). For interbranch 
collisions we need the relative energy e re i for a right and a left mover with 
relative momentum k. This is given by 

e re i (k) = -2E { p ] + 2hv F \k\. (24) 


For the interaction co, the principal value T-matrix is readily calculated. 
Notice that the inverse T-matrix is ultraviolet (UV) divergent, and we there¬ 
fore place a momentum cutoff A. The bare T-matrix takes the form T(T; k, k') 
r(T , )sgn(k)sgn(fc'), with 


r(e iel (k)) 


92 


1 


2nJivp lo § 


A—FI 

1*1 


(25) 


We now have two options. We can renormalise either the interaction cou¬ 
pling constant g 2 while letting the Fermi momentum take its physical value 
(kp^ = kp), or we may renormalise the bare Fermi momentum kp\ The 
renormalisation of the interaction coupling constant is straightforward. Ap¬ 
plying the renormalisation condition T(e re i(kp)) = g^\ we obtain 


1 

92 



1 

2nhvF 


log 


A — kp 

kp 


(26) 
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and then we have 


1 


r(e ie \(k)) 


27T hvp lo § 


k 

kp 


(27) 


The inverse T-matrix is now well-defined at almost all energies. It still retains 
an infrared (IR) divergence, which implies low-energy scattering is suppressed 
in this approximation. However, this is unimportant, since we will only be 
concerned with high-energy (of 0(2Ep)) processes. 

The other possibility, namely the renormalisation of the Fermi points, is 
quite easy to implement, too. Notice that for g 2 to be finite (and equal to 
g^), all we need is kp'* = A/2 + kp, with k F the renormalised (i.e. physical) 
Fermi momentum. In this case, the inverse T-matrix, Eq. (J25]), at energy 
2 Ep, is regular and equal to g 2 when A —* oo. The T-matrix then takes the 
form 


r(e re i(fc)) 


1 



2nk/Vp lo § 


1*1 


1 


(28) 


Written in the above form, the T-matrix appears to be ill-defined. However, 
if we are interested in relative momenta around the (bare) Fermi points, 

~ ~ / n\ 

k = A/2 + k, with k = 0(k F ), we have a constant T-matrix r(e re i(fc)) = g 2 
as we take the cutoff to infinity. For momenta far away from the (bare) Fermi 
points, we get a logarithmically suppressed effective interaction instead. 

All we need to do now is to obtain the renormalisation condition in terms 
of the phase shift at k — kp. In order for the renormalisation of the Fermi 
points to make sense, we need to assign the value of the phase shift at the 
physical value of kp to scattering at energy 2 Ep*. From the asymptotic anal¬ 
ysis of the scattering wave function we obtain the renormalisation condition 



—2Aui?tan 6k F , 


(29) 


which is exactly the same renormalisation prescription as for the quadratic 
dispersion, Eq. (|22]1 . showing that g 2 is unrenormalised when passing from 
the quadratic to the linear approximation. 

The renormalisation of the Fermi points, while leaving the interaction 
coupling constant be renormalised without infinities, naturally leads to Lut- 
tinger’s model. To see this, shift the momentum of right (left) movers as 
k —> k + A/2 (k — A/2), and rescale the cutoff as A/2 —y A. Then, the 
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momenta of right and left movers both belong to the interval (—A, A). The 
non-interacting Hamiltonian becomes 


h 0 — /iqn = hvp 


k =-A 


C kR^ C kR 


c kL^ C kL 


(30) 


where /zo = = fifvpA/A + Ep and N = IVr + Nl is the total number of 

particles, and where we have defined CkR (chif) as the annihilation operator for 
a right (left) mover with momentum k + A [k — A). The above Hamiltonian, 
Eq. (jS5P , exactly corresponds, after taking the cutoff A —> oo, to the non¬ 
interacting Luttinger Hamiltonian. This shows that, firstly, the renormalisa¬ 
tion of the Fermi points is a usual, implicit assumption in most treatments of 
Luttinger liquid phenomenology and, secondly, that the interaction coupling 
constant remaining ” unrenormalized” j23j (that is, renormalised without in¬ 
finities), is a consequence of the renormalisation of the Fermi points, and 
taking the limit A —> oo before allowing the coupling constant to change at 
all. 


5. Intrabranch processes 


The intrabranch coupling constant, typically denoted by g 4 23], is very 


simple to obtain as a function of the interbranch coupling constant, since it 
must obey a natural physical constrain, corresponding to the conservation 
of charge (particle number) [35?]. In our notation, it corresponds to setting 
g 4 = 2g\ . However, the understanding of the microscopic origin of g 4 is far 
from satisfactory, and we study this in detail in this section. 

It is customary to consider intrabranch interactions as two-particle scat¬ 
tering processes jH, 36| (see Fig. [2]), and then use a phenomenological cou¬ 
pling constant, namely g 4 , that simply adds to the Fermi velocity and renor¬ 
malises the bosonic excitation energies in the Luttinger Liquid [23[. There 
are, however, several problems associated with this collisional interpretation 
of intrabranch processes. Firstly, before linearising the single-particle disper¬ 
sion, one can define two-fermion scattering states in, say, the right-moving 
branch (we will consider right-moving fermions throughout this section un¬ 
less otherwise stated). Since the relevant scattering processes occur around 
the Fermi point, the energy at which one should renormalise the interactions 
should be twice the Fermi energy 2 Ep. However, the fermionic T-matrix 
(ergo the phase shift) tends to zero quadratically with relative momentum, 
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see Eq. (fSj) . This would mean that the renormalised Fermi velocity would 
not be changed to zero-th order in the momentum of excitations, contrary to 
well-known results and phenomenology [ 23 J. The second, even more severe 
problem, is that even if we wanted to renormalise the scattering properties 
in the linearised approximation to match the exact scattering properties (ob¬ 
tained using the full quadratic dispersion) around the Fermi point, we would 
simply not be able to do so. The fact, seemingly unnoticed so far, is that two 
particles (fermions or bosons) in the same branch of a linearised dispersion 
do not scatter. We begin by considering this issue in the following subsection. 

5.1. Absence of collisions 

Consider the Lippmann-Schwinger equation off the energy shell E ^ 0 in 
the relative coordinate, 

(k'\T(E)\k) = V(k,k') + f-V dq YlML( q \T(E)\k'), (31) 

where | k) is regarded here as a fermionic (antisymmetric) state, that is 
2 -1//2 (| k) — | —A:)) — > | k). As we go on-shell, E —)■ 0, we see that the 
principal value T-matrix must vanish linearly with E. Therefore, we de¬ 
fine T(E) = Et( 0) for E —> 0. As we take the limit E —> 0, we obtain the 
following integral equation 

1 

V (k,k') = -—J dqV(k,q){q\t(0)\kf). (32) 

The above equation can be written as an operator identity, which reads 

- V = Vt( 0). (33) 

The simplest way to solve Eq. (j33j) is by diagonalising V. We define a unitary 

operator f/@that diagonalises V, together with its associated orthogonal set 
of eigenfunctions |a*,) and their corresponding eigenvalues F(cKfc), such that 

K) = U\k). (34) 

Upon transformation of V and t(0) with U, Eq. (1331) reduces to 

- 27t V(a q )5(q - q') = V(a q i){a q >\t{ 0)|a 9 ). (35) 


2 Since V can be quite degenerate, the operator U is not necessarily unique. 
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Figure 2: Left: depiction and interaction diagram for intrabranch collisions. The red 
solid line corresponds to the Galilean (parabolic) dispersion, and the blue dashed line 
corresponds to the linearised approximation around the Fermi points. Right: depiction of 
the effect of the effective intrabranch interaction ( 34 ), to raise the excitation energies of 
quasi-particles (see text). 


For a strict finite-range potential V (which we can always consider as a good 
approximation), most of its eigenvalues will be degenerate and equal to zero. 
We must therefore deal with Eq. (j35]f with care. Let us define the set of 
values of q for which V(a q ) = 0, that we define as J 0 , that is 

h = U > 0 : V(a q ) = 0}, (36) 

while we denote its complementary set by Iy = 1 — h- If q G J 0 , then for 
each q' we must have either q' G J 0 or (a q '\t{())\a> q ) = 0. If q' G J 0 as well, 
then from Eq. (|35l) we see that (ccq/|t(0)|a 5 ) is arbitrary, which is a natural 
consequence of degeneracy. In this case, we shall fix its value to a constant 
co, i.e. 

(a g /|t(0)|a 9 ) =uj, q,q' G J 0 . (37) 

For q' G Iy , we are obviously forced to set (a q '\t(0)\a q ) = 0. We now consider 
the case q G Iy- From Eq. (l35ll we have that if q 7 ^ q' , then (cv^/|t(0)|«q) is 
arbitrary and we set it to the constant co as above if q' G Iq. The case q = q' 
yields (cvg|t( 0 )|cvg) = — 2tt. 

In the discussion above we have completely solved the Lippmann-Schwinger 
equation. We still need to calculate the scattering states. These have the 
form 

IV’) = !*:) + lim gW(E)T(E)\k). (38) 
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Inserting a resolution of the identity, we find the total scattering state as 


\i>) = \ k )+ 



dkf 

2 ^ 




It is very illustrative to choose u = 0 above. Using 
identity 1 in Eq. (l39]h and identifying 


- j i 7^( k '\ a q){ a ,\ k ) 

(39) 

the resolution of the 


1 _ Pio 



dq 

2tt 


otq)(ot, 


q\i 


(40) 


where Pj 0 is the projector onto Jq, we End that 


= pm. 


(41) 


The above result has a clear interpretation: the scattering state of the system 
associated with the incident wave | k) is nothing but the projection of the 
incident wave onto the zero potential eigenstates. An immediate consequence 
of Eq. (l4Tf is that our system (or any other flat-banded system) only supports 
scattering states if the potential has at least one vanishing eigenvalue. 

The above analysis proves our claim that intrabranch interactions can 
lead to either (i) energy shifts if the potential has no zero eigenvalues or (ii) 
a superposition of non-interacting waves that avoid potential overlap if the 
interaction has vanishing eigenvalues E|. 


5.2. Fermion-fermion intrabranch interactions 

It is clear from the discussion in the previous subsection that the re¬ 
processes do not correspond to particle-particle intrabranch scattering. Lut- 
tinger liquid theory, on the other hand, predicts that particle-hole excita¬ 
tion energies are shifted due to intrabranch interactions. We will see in the 
following that studying the energy shifts of particle-hole excitations in the 
chiral (right- and left-moving) branches separately allows for some sort of 
renormalisation of the coupling constant in the weak inter -branch inter¬ 
action regime. We will show, however, that microscopic fermion-fermion 
intrabranch interactions do not correspond to the intrabranch interactions 


3 These correspond to the ’’strange solutions” of one-dimensional field theories found 
by Mattis and Sutherland 421. 
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appearing in Luttinger liquid phenomenology. These will be considered in 
the next subsection. 

We begin by defining the chiral right-moving problem, with single-branch 
Hamiltonian H r given by 


H R = Y,<k)c] c c k + V BR , ( 42 ) 

k> 0 


where 

Vrr = ^L V ( Q) c k+A’- q c k'CkO(k + q)9{k' - q ). (43) 

k,k'>0,q 

Hamiltonian (142|) conserves total momentum. Since the Fermi sea \Nr) for 
Nr right-moving fermions, given by 

tVfl>= n 4|0>. (44) 

0<k<kp 

is the unique state with lowest total momentum, it is an eigenstate of the 
interacting Hamiltonian Hr. For quite general repulsive interactions we can 
also assume it is the ground state of the chiral system. In order to obtain 
particle-hole excitation energies, it is convenient to rewrite the interaction 
(l43li in the following way (throwing away constant energy shifts) 

VrK=J^ V (q)Pq p q ^) C l+qCk+q (45) 

g>0 kq> 0 

Above, we have defined the displacement operators p q R and p+ R as 

PqR = C k-q C k 

k>q 

PqR — C k+q C k = P~q 

k> 0 

Since the first term on the r.h.s. of Eq (145|) is bosonizable, we shall first 
deal with the second term, which we denote by V. In particular, we are 
interested in the difference in energy between particle-hole excitations and 
the Fermi sea due to this interaction term. We will denote particle-hole states 
by | k + Q-,k), with 

\k + Q;k) = cl +Q Cj.\N R ), (48) 


( 46 ) 

( 47 ) 
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where obviously k < kp and Q > 0. The action of V on a particle-hole state, 
throwing away a trivial constant term so that it gives the excitation energies, 
is given by 


V |k + Q ; k) 


1 

L 


k+Q 


V(k - k) - J2 V(k + Q - k) 


k =0 


k =0 


| k + Q-, k), 


(49) 


which means particle-hole states are eigenvectors of V. Since we are interested 
in Luttinger liquid phenomenology, we will consider only excitation energies 
of 0(Q). To this order, we can obtain these analytically by expanding V[k + 
Q — k) as 


VCk + Q-k ) 


V(k-k) + 


dV 

dn 


Q + 0(Q 2 ). 

K=k—k 


(50) 


We obtain 


V\k + Q ; k) 



V(Q) — V(k) Q 


1 

27T 


' 'k+Q 


dkV (k — k ) 


In order to obtain the leading order in Q above, we set V ( Q ) 
and finally arrive at 


• | ~k + Q-,k)+0(Q 2 ). 

(51) 

V(0 ) + 0{Q 2 ) 


V| k + Q ; k) 


1 

2tt 


V (k)Q + 0{Q 2 ) 


| k + Q] k). 


(52) 


Above, we see that the excitation energy to linear order in Q due to the part 
V of the intrabranch interaction is given by —V{k)Q/2Ti. 

We now consider the first term on the r.h.s. of Eq. (l43lh which we 
will denote by Vb , he. the part of the interaction that is straightforwardly 
bosonizable. For sufficiently small q (q -C kp), the usual commutation rela¬ 


tions for the displacement operators 23] are valid, 


\PqRi P^'R\ 



(53) 


We define bosonic operators a q R = (2n / Lq) 1 ^ 2 p q R (and equivalently for a] y7? ), 
and the term Vb is rewritten as 

Vb = 2 ^ qV ^) a lR a iR- ( 54 ) 

g>0 
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The eigenstates of Vb including only single particle-hole pairs are given by 



and their respective eigenvalues are given by qV( 0) + 0(q 2 ). 

In order to establish the excitation energies solely in terms of a quadratic 
bosonic Hamiltonian, we need to make two approximations at this point. 
One of them is the linearisation of the single-particle dispersion, so that the 
excitation kinetic energy of a particle-hole state \k + Q\ k), which is given by 
h 2 kQ/m + 0(Q 2 ), is approximated by hvpQ/m for all k. The second one 
corresponds to approximating the energy shift due to V as 

- h v(k)Q “ ~h v(kF ) Q ' (56) 

for all k, a reasonable approximation for small (). In that approximation, 
a \ R \Nn) are also eigenstates of V, and we can finally write down the intra¬ 
branch interaction, approximately, as 

Vrr = ^Y.1 iC(°) - C(M] d\ R a, R . (57) 

g>0 

The above intrabranch interaction would give a (weak-coupling) value of 
g 4 = V(0) — V(kp)- This is incorrect, since in the weak-coupling limit the 
result should be g± = V(0) — V(2kp). What has gone wrong here is that the 
intrabranch interaction in the effective Luttinger’s model does not correspond 
to an intrabranch process in the original fermionic system. We will see how 
this is the case in the following subsection. 

5.3. Luttinger liquid intrabranch interaction from the original interbranch 
process 

We consider now the weak-coupling limit of the intrabranch interaction. 
The relevant state to be considered is a particle-hole excitation from the 
static Fermi sea in the original fermionic system, 

I k + Q)= 4 +q c s |F), (58) 
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where 0 < k < k F and 0 < Q <C k F , while 


kp 

i f )= n 4io). 

k=—k p 

We denote the relevant part of the interaction by V', which is given by 
V = 7E b(°) -V(k-k-Q) 1 4c4 +0 c*. 


( 59 ) 


(60) 


All other interaction processes give contributions of 0(Q 2 ) or higher. We 
obtain 


V’\k + Q) = 


2 tt 


dk 


—kp 


V(0)-V{k-k-Q) I k + Q,k). (61) 


Subtracting the energy shift corresponding to the Fermi sea, i.e. that ob¬ 
tained by setting Q = 0 above, and after using k ~ k F , we obtain the energy 
shift g^Q/2-n to O(Q), with 


4 1 ’ = v(o) - v (2 k F ), 


(62) 


which is the correct answer we were looking for. 


6. Final Flamiltonian and Luttinger parameters 


With the notations used in the previous sections, we can write down the 
effective Luttinger liquid Hamiltonian of the system in bosonized form, i.e. 
in terms of the boson operators for right (left) movers a qF ( a q L )• This is 
given by 


H = hv ^ q 

q> o 


a \il a qR + a li a qL 


+ 


92_ 
2 vr 


E" 

q> 0 


a qR a qL + a qLRqR 


(63) 


where v = v F + g 4 /27rh, and where #2 = 2(7^, see Eq. ()29]h The last 
statement might not be so easy to grasp, so we prove it as follows, using the 
notation of ref. (35|. The original interbranch interaction is given by 


v<2) = jz E [ y («) - v & -*-«)] 

k>0,k'<0,q^0 

X C kR C k'+q,R c k’,L C k'-q,L + (R L). (64) 
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To see how the p-wave interaction arises due to the exchange above, we 
define Q = [k — k')/2 + q and Q' — (k — k')/2 , in which case q = Q — Q' and 
k' — k — q = —(Q + Q 1 )- Therefore, we have V(q) — V{k' — k — q) = 2V P (Q,Q'). 
Replacing now the p-wave potential by the ’’constant” interaction, we obtain 

Vp(Q, Q') = g 2 sgn —— + q^j sgn ■ (65) 


Then for the first term in Eq. (IR-fl) we have (k — k')/2^kp>0 and, for small 
q (forward scattering), we have + q > 0 , and therefore V P (Q,Q') ~ g 2 , 
which leads to the ^-interaction term in Eq. (1631) . 

Hamiltonian f|63)l corresponds to the notation in ref. j23|, and the Lut- 
tinger parameter /C and the speed of excitations u are given by 


K, 


^ 1 9 4 _ 9 2 

2'kKvf 2 ttHvf 

1 _L 94 I 92 

27 tHvf 2'kKvf 


1/2 


U = Vf 



g 2 

2nhvF 


21 V 2 


( 66 ) 

(67) 


In the weak-coupling limit, we can expand the Luttinger parameter and the 
speed of excitations as follows 

IC = 1 - 2 ^ + O(sl), ( 68 ) 

“ = VF 0 + 2 ^) + ° fe2) ' (69) 

Our simple renormalisation prescription gives the first order results above 
exactly. We have to notice that the first order results in g 2 (or equivalently 
( 74 ), obtained from our renormalisation prescription are not perturbative, 
since they represent the exact two-body phase-shifts of the system and not 
their first Born approximation oc R(0) — V(2kp). 


7 . Lieb-Liniger and Cheon-Shigehara models 

We apply now the results of the previous sections to a particular example 
of an exactly solvable model, namely the (bosonic) Lieb-Liniger gas jl3|, or 
its fermionic dual, the Cheon-Shigehara model [38| . 
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For the sake of clarity, and to avoid confusion with other notations in 
the literature, we write down the Lieb-Liniger Hamiltonian in first 

quantised form 

H Li eb = + 9 6 ( Xi ~ x ^- ( 70 ) 

i i<j 


The dual fermion-fermion interaction in Cheon-Shigehara’s model, however, 
has a formidable form in the position representation. Its momentum rep¬ 
resentation, on the other hand, is very simple, and does not appear to be 


/ery 

widely known. This is given in ref. 43|, and has the form 


Vcs(k,k') = g F kk', 


(71) 


which corresponds to a p-wave interaction (see Eq. ([ 6 ])). From the interaction 
in Eq. (ITT]) . we obtain the following phase shift [43] 


, Q rng F 

tallies = — ypk. 


(72) 


The dual interaction in Lieb-Liniger’s model is obtained by equating the #cs 
to the Lieb-Liniger phase shift 13], obtaining 


4 fh 2 


9 f — -I — 

g \m 


(73) 


The above discussion immediately yields the value of $2 (and gq) in terms of 
the Lieb-Liniger coupling constant g, 

92 = ~~(hv F ) 2 . (74) 

9 

Defining now the dimensionless Lieb-Liniger parameter 7 , 

_ 7T mg _ -Kg 
7 k F h 2 hv F 

we obtain the desired first order result in the Tonks-Girardeau limit 

/C = 1 + - + e>( 7 - 2 ) (76) 

7 

n = l-- + 0( 7 - 2 ), (77) 

7 
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in agreement with the known result (see, e.g. jUj). 

It is worth noting that in this model only , the first Born approximation 
to $2 is identical to the result using the phase shift. This is due to the fact 
that Lieb-Liniger’s model is a zero-range S-matrix theory, i.e. the interaction 
has zero range and the model is solvable via the Bethe ansatz. In collision- 
theoretic language, the interaction potential in this model has already been 
renormalised to give the correct phase shift. To see this, we note that the 
following interaction has the same action on fermions as the Cheon-Shigehara 
interaction CD, 

V(q) = ~q 2 , (78) 

which implies V (0) — V(2k F ) = —(8 /g)(hvp) 2 , which is identical to the result 
obtained using the phase-shift, Eq. (l74l) . 

We illustrate now the fact that the first Born approximation to g 2 being 
identical to the non-perturbative result only occurs in Lieb-Liniger’s model, 
where the interaction has been renormalised already for that purpose. We 
consider a full-blown two-particle interaction of the following form in the 
position representation 

W(x) = g 0 e~ x ^ : (79) 

where x is the relative coordinate, go is the strength of the interaction, and 
A > 0 controls its range. For two fermions, the phase shift is readily calcu¬ 
lated from the relation 

/ OO 

dyW(y)sm(ky)i> k (y), (80) 


where the fermionic scattering state ifk(x) satisfies the Lippmann-Schwinger 
equation 


V’fc(z) 


sin (kx) + 


m 

2 ¥k 


dy sin( 


x 


y\)W(y)ij}k(y)- 


(81) 


The results for the phase shift (which give the value of gf) are shown in Fig. 
[3], where we have set A/fcp = 1. There, we observe that the weak-coupling 
result is only valid in a small window of weak interaction strength go/Ep. 
However, the gas phase of the system can be weakly-coupled for strong at¬ 
tractive interactions, to the left of the shape resonance in Fig. [3j In fact, 
the strongly-attractive side of Lieb-Liniger and Cheon-Shigehara models (the 
super-Tonks regime), corresponds to this region in the exponential potential’s 
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9 /2 f 

Figure 3: The tangent of the phase shift (solid red line) at relative momentum k = Rf 
for two fermions interacting via an exponential potential, Eq. CZSD, compared to its first 
Born approximation (dashed blue line), for X/kF = 1. The vertical red line denotes the 
position of a shape resonance. 


strength, while the weakly attractive regime in Cheon-Shigehara model cor¬ 
responds to the weakly attractive regime of the exponential potential, and 
so does the strongly repulsive (Tonks regime) Lieb-Liniger model. We also 
note that there is no such thing as repulsive Cheon-Shigehara model. When 
gF > 0 (g 2 > 0) we still have a strongly-attractive interaction due to renor¬ 
malisation j43j. Cheon-Shigehara model is everywhere attractive: strongly 
attractive (there are bound states) for g F > 0 (g 2 > 0) and weakly attractive 
(there are no bound states) for g F < 0 (g 2 < 0). 

8. Conclusions 

In this article we have studied how the phenomenological constants of Lut- 
tinger liquid theories relate to the microscopic interactions of spin-polarised 
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one-dimensional Fermi gases. We have shown that interbranch interactions 
rule both phenomenological inter- and intrabranch interactions in Luttinger’s 
model. This has also been shown to arise in a completely natural way from 
the renormalisation of the Fermi points in Tomonaga’s model, upon an ap¬ 
propriate redefinition of the wave numbers for right- and left-movers. As 
an important side result, our renormalisation programme explains, in very 
simple terms, the well-known fact of unrenormalisabilty of the interbranch 
coupling constant. 

Our theoretical framework can be generalised to Luttinger liquids with 
spin, chiral Luttinger liquids - where the results of section 15.21 are directly 
applicable - and helical Luttinger liquids j45[ . Many of the concepts studied 
here can be relevant for the study of non-linear Luttinger liquids [44|, i.e. 
when the non-linearity in the kinetic energy dispersion becomes important. 
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